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Chapter 1 
Reaction mechanisms of pair transfer 



Ian J. Thompson 

Lawrence Livermore National Laboratory, 
PO Box 808, L-4U, Livermore, CA 94551, USA 

The mechanisms of nuclear transfer reactions are described for the 
transfer of two nucleons from one nucleus to another. Two-nucleon 
overlap functions are defined in various coordinate systems, and their 
transformation coefficients given between coordinate systems. Post and 
prior couplings are defined for sequential transfer mechanisms, and it is 
demonstrated that the combination of 'prior-post' couplings avoids non- 
orthogonality terms, but does not avoid couplings that do not have good 
zero-range approximations. The simultaneous and sequential mecha- 
nisms are demonstrated for the 124 Sn(p,t) 122 Sn reaction at 25 MeV us- 
ing shell-model overlap functions. The interference between the various 
simultaneous and sequential amplitudes is shown. 



1. Introduction 

Much of the evidence for nucleonic pairing in nuclei comes from energy ex- 
pectation values, but important further information comes from the transfer 
of pairs of nucleons to or from another nucleus of known structure. In this 
regard, a more fundamental understating of nuclear reactions has been, 
and will continue to be (especially in the FRIB era), crucial to the nuclear 
physics community. This chapter focuses on the theory, calculation and 
model results for the reactions mechanisms of pair transfer. 

Here we consider the reaction mechanisms for pair transfer between two 
nuclei, namely reactions that we can describe as A(B+2, B) A+2. Here, the 
two nucleons may be two neutrons, two protons, or a proton and a neutron, 
and are transferred from core B to core A. The nucleons may transfer either 
in one simultaneous step, or one after the other sequentially. If a distinguish- 
able proton and a neutron are transferred, then both proton-then-neutron 
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and neutron-then-proton routes need to be considered. Furthermore, these 
sequential and simultaneous routes contribute amplitudes that all add to- 
gether coherently. This feature enables us to probe the nature of coherent 
two-nucleon superpositions in nuclei. Conversely, these superpositions, cou- 
pling orders and phase conventions have all to be defined consistently in a 
good calculation. 

Subsequent sections will therefore consider the definition of two- 
nucleon overlap functions, their coordinate transformations, the definition 
of transfer matrix elements along with zero-range approximations and non- 
orthogonality corrections. Finally, some results are shown to illustrate the 
coherence effects in the reaction mechanisms of pair transfers. 

In the last 50 years, a significant number of papers have been presented 
in which absolute differential cross sections have been calculated, and com- 
pared with experimental results. 1 " 19 Traditionally, for example inj^ the 
theory predictions have fallen well below the experimental data. This ratio 
has been called the 'unhappiness factor'] 20 ! 21 ! and has sometimes been more 
than 100. Most previous calculations modeled the transfer of a dineutron 
as a single cluster. And only from CharltorP were sequential transfer con- 
tributions considered. We find that modern calculations, such as as those 
of Potel et alp^ are in considerably better agreement with experiment. 

2. Bound states and vertex functions 

The general theory of nucleon pair bound states defines the overlap function 
<pj(r,p) — (&a(I)\&a+2(J)} in terms of the Jacobi coordinates r between 
the two nucleons, and p between their center of mass (cm) and the core A. 
The core spin is I and the spin of the ^4+2 composite state is J. When 
intrinsic spins S\,S2 are also included in a particular coupling order such as 
\{L, (£, (sis 2 )S)j} J12, 1', J), we have the partial wave expansion 



<t>i{v,p)= Y, "W&OCC Yt$)Y?{p) -u 12 (r,p)(J 12 M 12 I^\JM) 



The radial wave function u\ 2 (r, p) can be given either as a cluster product 
of single-particle wave functions u\ 2 {r,p) = <&i,{r)4>t(p), input directly as a 
two-dimensional distribution e.g. from a Faddeev bound-state calculation, 
or calculated from the correlated sum of products of single-particle states 
with independent coordinates. These two-nucleon wave functions will in 
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general be the eigenstates of a three-body bound state Schrbdinger equation 

[T r + T p + Via + V 2A + V 12 - e]0/(r, p) = , (2) 

where the Via are the potentials between nucleon % and the core, and V\ 2 
is the pairing interaction between the two nucleons. 

Such two-particle states that come from shell-model calculations^! or 
from Sturmian-basis calculations^! are usually described by means of the 
|ri,r2) coordinates. This describes a pair state as 

Vi2(ri,r 2 ) = \(li(i),s 1 )j 1 (i),(£ 2 (i),s 2 )j2(iy, J12) (3) 

i 

The coefficients Ci for correlated basis states % and the single-particle wave 
functions (f£ S j (r) contain all the physics information about the bound state 
needed to do a transfer calculation. Shell model codes^! can produce the 
coefficients Cj needed here in terms of previously calculated eigenstates of 
the A and the A + 2 systems. These states are then transformed into the 
centre-of-mass coordinates |r, p) of Eq. (fTJ) using Yi = Xir + yip. For equal 
mass particles, x% = x 2 = 1, and yi = —y 2 = \. 

The vertex functions of these bounds states are defined to be these 
bound state wave functions 0/(r, p) multiplied by the potentials which 
have zero effects after the transfer step is performed and all exit channel 
nuclei have completely separated. These potentials are therefore the sum V 
of the binding potentials V^ sp = Ve S j(ri), namely V = V^ + V 2 S q. (These are 
the individual potentials that should appear in the bound-state equation 
[T r + Ve S j(r) — e](fi£ S j(r) = 0.) The vertex function does not include the 
nucleon- nucleon pair interaction Vi2(i"i— r 2 ), since this potential produces 
binding effects in both the initial and final bound states. We denote by 
V4>j{y, p) the vertex function after transformation into Jacobi coordinates 
by the same method used to transform the wave function itself. 

3. Post and prior coupling forms of transfer matrix elements 

We now consider the Hamiltonian W for the whole system of A+B+2 nu- 
cleons and described by system wave function ^ for the complete transfer 
reaction A(B+2, B)A+2. Let the various A+2 and B+2 bound states be 
denoted by $i for indices i. Then we may expand ^ in terms of the $i with 
some coefficients ^(R^) depending on the two-body separation vectors R^. 
This gives the channel expansion ^ = ^ ^(R^)^. 
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The model Schrodinger's equation [H — E]^> = when projected sepa- 
rately onto the different basis states $j , yields the set of equations 

[Ej - Hj] tpj(Rj) = ]T (^\n - E\$i) ^(Ri), (4) 

which couple together the unknown wave functions ^(Rj). The channel 
Hamiltonians are defined by the diagonal Hj — Ej = ($>j\W — E\$>j) such 
that the Ej are the asymptotic kinetic energies in channel j. 

The off-diagonal matrix element (<&j\H — E\$i) has two different forms, 
depending on whether we expand 

H — E = Hj — Ej + Vj (the 'post' form) 
= Hi — Ei + Vi (the 'prior' form). 

The name (post or prior) is determined by whether it is the initial or final 
channel whose Hamiltonian is used. The above Eq. (j4]), as written, has i 
as the initial channel and j as the final channel for the indicated coupling. 
Thus 

(QjlH - E\$i) = ^ ost + [Hj ~ E^Kji (post) (5) 
or = V}f or + Kji [Hi - Ei] (prior) , 

where 

Vfr^^mm, ^ rior = <^im>, Kji^i^i). (6) 

The overlap function Kji = ($j|$i) in Eqs. (15161) arises from the non- 
orthogonality between the basis states 3>j and $j if these are in different 
mass partitions. The Kji are non-local operators that go to zero asymp- 
totically {Ri or Rj —¥ oo). (Within the same partition, the $i would be 
inelastic states, and would form an orthogonal set.) 

The first-order DWBA matrix element use entrance ipi and exit if)j chan- 
nel wave functions satisfying [Hi — Ei]ipi = and [Hj — Ej)ipj = respec- 
tively. Its matrix element is 

Tj? = {^- ) $ j \H-m4 + )) . (7) 

The prior form of this is 

T (prior) = { ^-)^ H . _ E . + Vi\^4 +) ) 

= (^-^j\Vi[^i4 +) ) + {^QjlHHi - Ei}4 +) ) 

= (^\Vf OI \4 + ^) . (8) 
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Similarly, the equivalent post form is 

H^ _) ^-lM +) > + o 

= ^ j - ) \vfr t \^ ) ) ■ (9) 

Thus the non-orthogonality term disappears in first-order DWBA. Post 
and prior first-order DWBA matrix elements can be made to exactly agree 
numerically, if sufficient care is taken to ensure convergence of the non-local 
form factors. 

Let a second-order DWBA matrix element use entrance channel i, exit 
channel k, and some intermediate channel j, as i — > j — > k. The propaga- 
tion in the intermediate channel may be described in terms of the Greens 
function Gj, or equivalently within an iterated coupled-channels set. The 
two-step DWBA matrix element is 

= (4~ } \ {^\H-E\^ 3 )G 3 ^ 3 \H-E\^) . (10) 

Now there are four matrix elements that may be calculated, according to 
the first and the second Hamiltonian form chosen: post-post, post-prior, 
prior-post, and prior-prior. The terms prior and post for each step are used 
to refer to the initial or final channels of that step, not the overall incoming 
or outgoing channels. In 'prior-post', the prior refers to the first step, and 
the post refers to the second step. 

The post-post form of this, for example, is 

T (post, P os t) = ^ klHk _ Ek+Vk ^ j)Gj ^ jlHj _ Ej+V . m ^(+)). 

Here the [Hk — Ek] can operate on the final tpk to give zero, but little can 
simplify the [Hj — Ej] since [Hj — Ej]Gj 7^ always. Thus 

T ( V ost, P ost) = ^(-) |y post G ,. F po S t| V) (+) ) + 

(4~ } \K; stG ^ - ^ (+) >- (12) 

This second term is called a 'non-orthogonality term' since it involves the 
bound-state non-orthogonality overlaps Kji = which is significant 

when Ri and Rj are both within the range of the bound states. 

Similar analyses for post-prior and prior-prior two-step DWBA expres- 
sion also have non-orthogonality terms in the final form. The prior-post 
form, however, is 

T ( pri or, P ost) = mHk _ Ek+Vk ^. )G .^.\ H ._ Ei+V ^. ) ^ + )).(13) 
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Here the [Hi — Ej] can also operate on the initial ipi to give zero, as well as 
[Hk — Ek] on the ipki so we have the simplest form 



The non-orthogonality terms can thus be made to disappear in second-order 
DWBA if the first and second steps use the prior and post interactions 
respectively. If the non-orthogonality terms are included as necessary, the 
results should be the same whatever post or prior forms are used. 

In third and higher-order transfer calculations, some non-orthogonality 
terms will always be present, but most pair transfer mechanisms can be 
well modeled as two-step processes. 

4. Two-nucleon transfer interaction 

We now consider the specific transfer matrix element V^ rlor = ($j|Vi|<I>i). 
Given an expression for this prior form, we may calculate the post inter- 
action easily as V^° st = (V? llor )^ . Take $j to refer to the bound states of 
nucleus A+2 outside core A, and $i analogously for nucleus B. 

The transfer interaction has therefore the non-local matrix element 



V J4 (R 3 , R0 = (r, p A )\Vll + V% + U AB - (r, p B )). (15) 



As is usual in transfer operators, there are three kinds of potentials appear- 
ing here. First there are the binding potentials ^^(ri^ + V^^s). Since 
these binding potentials always appear while multiplied by their bound 
state wave functions, we need only store and use the vertex functions de- 
fined in section [21 Secondly, there is the 'core-core' potential Uab(R-ab) 
between the core nuclei A and B. Finally is subtracted an optical potential. 
In this prior form we subtract the optical potential in the initial channel, 
E/j(Rj). The difference U ab — Ui of the two optical potentials is called the 
remnant term, and is sometimes taken to be small. 

The integrals in Eq. (|15j) include integrating over the two-nucleon sepa- 
ration r as well as over their cm distance p A from the core A. The r coordi- 
nate appears in both the initial and final states, and so is not labeled by A 
or B. This has the important consequence that neither the distance nor the 
angle of the r coordinate is changed in the transfer. Neither, therefore, is 
their relative angular momentum £, and, for similar reasons, nor their spin 
couplings S and total angular momentum j. The two neutron transfer can 
hence be viewed as the transfer of a 'structured particle' {r, (£, (siS2)S)j}, 
and then becomes similar to the more familiar single-particle transfers. This 



T, 



-t(prior, post) 
hi 



(14) 
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means that when we also afterwards integrate over the coordinate p A , we 
can use the standard procedures already developed for one-particle transfer 
interactions. 

5. Coordinate transformations 

The transfer mechanism requires the pair wave function to be expressed in 
the form of Eq. ([1]), so independent-particle forms of Eq. (O have to be 
transformed in their coordinates as 

Viagra) = J> £ \L, (£, ( Sl s 2 )S)j; J 12 T) 0* ( ££(r,p) . (16) 

u LISj 

A particular basis state % in the (r, p) coordinates is 
<t>L{M)j( r >p) = (L,(£,(s 1 s 2 )S)r,J 12 \(l 1 ({),s 1 )j 1 (i),(£ 2 (i),s 2 )h(i); J 12 ) 
x([Y L (r)Y e (p)} x | [w lSUl (ri)w !S232 (r 2 )] Ji3r ) (17) 

where (suppressing the i indices for clarity), and including an isospin T to 
define the antisymmetrization, 

{L, (£, (sxs 2 )S)j; Jx 2 T\(£x, s\)j%, (£ 2 , s 2 )j 2 ; J\ 2 T) = ^ XSj%j 2 

A 

^1 ^2 A \ ( y+s+t 

Sl s 2 S ' = = jXW(L£J 12 S;Xj)(-l) e+L -\ (18) 



V2(T + Se 1 ,e 2 S jl ,j 2 ) 



The radial overlap integral can be derived by means of harmonic-oscillator 
expansions with the Bayman-Kallio expansion^ or using the Moshinsky 
solid-harmonic expansion.^ This last method gives 

Kkt A (r,p) = ([Y L (r)Y e (p)]x I IwMwA^t) (19) 

1 1 



n±n-2 
X 



2ni J \ 2n 2 
X)q&,(r,p) (2Q+1) £ x £ 2 £^ x £^ 2 L£ 



li-m n 2 A x \ / C-n 2 n x A 2 \ (Ax L Q\ f A 2 £ Q 

o o o J I o o o j o oo) 00 



*E 

AiA 2 

x(-l) £l+£2+L+A2 (2A 1 + 1)(2A 2 + 1)W(A X LA 2 £; QX) 
£ x —n x n 2 Ai 

x I m £ 2 -n 2 A 2 1 , (20) 
£i £ 2 A 
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where y^j * s a binomial coefficient. The kernel function q^ 2 (r, p) is the 

Legendre expansion of the product of the two radial wave functions in terms 
of tp, the cosine of the angle between r and p: 

I I \^2 + l 



i r +1 



n r 2 



P Q {u)du (21) 



6. Zero-range and other approximations 



The coupling potentials Vj^R^Rj) of Eq. (|T5|) are non-local, in the sense 
that in general the initial and final radii, Rj and Ri, will be different. They 
will not only have different magnitudes, but also different directions. In 
the early days of transfer modeling, the calculations only became practical 
if a zero-range approximation could be found, in which the coupling was 
restricted to Rj = oRj for some constant a (which need not be unity). 

When the projectile is a light ion such as 3 H, 3 He or 4 He for nucleus 
B + 2, then the binding potential sum V? B + V^ B will have short range. We 
may therefore consider approximating the vertex function 

MS + ^Ms^Pb) ~ D a 5(p B )^ n (r) (22) 

for some nucleon-nucleon wave function (j> nn (r) that we are free to choose. 
This a zero-range approximation. Note that it is only p B which needs to 
have zero range, not r. The constant D$ is called the zero-range constant. 

If, furthermore, we can neglect the remnant term Uab — U{, then the 
transfer coupling of Eq. (|15l) can be simplified as 

v,.(R,.R,> = (4>%(t, Pa )\v& + v£\4>%(t,p b )) 

= (<t>%(r,p A )\D 6(p B )<P% n (r)) 
= D (<^ n (r)|<^(r,pJ) 6(p B ) 

= D Q (^ n (T)\^(r )PA )) 6{/3 (R 3 ^ R *))> (23) 



since 



Rj r— Ri = Pb/P for P= 2{ A A \ B+2 \ - (24) 

J A+2 yB/H M (A+2)(B+2) v ' 

That is, we arrive at a 'form factor' (^ n (r)|<£j (r, — Rj)) that is local in 
Rj = -^p^YLi = — p A because of the delta function 8{p B ). To find the form 
factor, we need to determine the average nucleon-nucleon relative wave 
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Table 1. Two- neutron overlap 
function for ( 122 Sn| 124 Sn) 

lg 2 /2 0.62944 

2d 2 /2 0.59927 

2dl /2 0.71913 

3s 2 /2 0.51892 

lfi. 2 1/2 -1.24399 



function <p^ n (r) in the light ion, and project the heavy- nucleus two-body 
wave function <pj^(r, p A ) onto this relative motion. This gives a function 
only of the distance p A = Rj and the angles. The kinematics for this zero- 
range approximation are identical to those for the one-body transfer of a 
mass-2 cluster from core B to core A. A local-energy approximation may 
be used to improve the treatment of the finite range of the vertex function, 
just as for one-body transfers. 

This is a further example the conclusion stated at the end of section 
[21 namely that transfer reactions only probe in the unknown nucleus those 
components of nn relative motion that already exist in the known nucleus. 
Since the known light nuclei 3 H, 3 He and 4 He have predominantly s-wave 
relative motion between the two transferred nucleons, our transfer reactions 
will only probe pairing states of s-wave relative motion in the target. The 
magnitude of the transfer cross section will be proportional to the form 
factor overlap (<^ n (r)|0^(r, p A )). 

Zero-range approximations can be also used for some of the sequential 
steps involving these light nuclei, but not for all of them if we are using 
'prior-post' couplings to avoid non-orthogonality corrections. For stripping 
reactions such as (t,p), the first prior (t,d) step has no good zero-range 
approximation, and for pickup reactions such as (p,t), the second post (d,t) 
step must be treated in full finite range for the same reason. 

7. Results 

In this short paper we will focus on the reaction mechanisms for the pair 
transfer 124 Sn(p,t) 122 Sn at 25 MeV, using the overlap function shown in Ta- 
bled] we find by overlapping the shell- model wave functions for the ground 
states of 122 Sn and 124 Sn. The structure results for 124 Sn and 122 Sn were 
obtained in the model space of (0*77/2 , W5/2, l<^3/2> ^ s i/2i O/in/2) for neu- 
trons with the code NuShellP^ The model-space two-body matrix elements 
arc those used in RefsPmpj They 

were obtained starting with a G matrix 
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Fig. 1. Simultaneous (short dash), sequential (dot-dash) and simultancous+sequential 
(solid line) cross sections for the reaction 124 Sn(p,t) 122 Sn at 25 MeV, in comparison with 
the experimental data of Guazzoni et a/EIl 



derived from the CD-BonrPHnucleon-nucleon interaction. The harmonic os- 
cillator basis was employed for the radial wave functions with an oscillator 
energy hw = 7.87 MeV. The effective interaction for the above shell-model 
space is obtained from the Q-box method and includes all non-folded di- 
agrams through third-order in the interaction G and sums up the folded 
diagrams to infinite order P2M1 The single-particle energies were adjusted 
to reproduce the observed states in 131 Sn. 

The inputs to the reaction code are the two-nucleon spectroscopic 
amplitudes (TNA) of Table [TJ A center of mass correction^ equal to 
[A/(A-2)] 2n+e for the TNA has been applied, where A = 124. Our sign 
convention is that the radial wave functions are positive at the origin. The 
sequential process was calculated by a single intermediate state for each 
of these orbits connected by a product of one-nucleon spectroscopic ampli- 
tudes that are equal to the center-of-mass corrected TNA multiplied by \/2 
that takes into account the normalization of the two-particle amplitude. 
Future calculations should also take into account the TNA obtained from 
the mixing of neutron pairs for orbitals outside of the model space. 
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Scattering angle (degrees) 

Fig. 2. Sequential cross sections for all possible combinations of post and prior for the 
two steps. 

We use the triton potential of Li the deuteron potential of DaehnickJ^ 
and the proton potential of Chapel Hill 89 All the two- neutron wave 
functions are constructed within the half-separation-energy prescription. 
For a triton wave function we use the pure s 2 configuration found by the 
product of eigenstates at the half-separation energy (4.24 MeV) in a Woods- 
Saxon potential with V = 77.83 MeV, R = 0.95 fm, and a = 0.65 fm (the 
results are not sensitive to these values). The Sn wave functions shown 
in Table [T]are found at the half-separation energy (7.219 MeV) in a WS 
potential with r = 1.17 fm, and a = 0.75 fm that has the fixed spin-orbit 
component V so = 6.2 MeV, r — 1.01 fm, and a = 0.75 fm. 

The complete cross section prediction is shown in Fig.[TJ compared with 
the experimental data of Guazzoni et a/PSl Now we see that, with the shell- 
model overlaps and proper finite-range and sequential contributions, the 
unhappiness factors are much closer to unity. A better agreement between 
theory and experiment has already been published,^ but in the present 
calculations there are still questions about the angular oscillations which 
are in not so good agreement with experiment. Note that Guazzoni et al^^ 
took the better agreement of the simultaneous transfer curve (dashed line) 
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Fig. 3. Simultaneous (short dash), sequential (dot-dash) and simultancous+sequential 
(solid line) amplitudes at zero degrees for the reaction 124 Sn(p,t) 122 Sn at 25 MeV. The 
short lines show the individual contributions from the wave function components of Table 
[T] and the longer lines with symbols are their coherent sums. 

to indicate small effects for sequential transfers, but this is not correct since 
we do know that sequential transfers occur, and can calculate them with 
good accuracy in this model (dot-dashed line). 

To see the importance of the non-orthogonality terms, and hence of 
choosing 'prior-post' couplings if non-orthogonality terms are to be avoided, 
Fig. [2] plots the different sequential cross sections for all possible combina- 
tions of post and prior for the two steps. The prior-post solid curve is the 
dot-dashed curve in Fig. [TJ The other curves are all different from this 
one, and cannot be simply added to the simultaneous amplitude to get the 
correct result. This also implies that no complete calculation with only 
zero-range couplings is possible. 

Finally, it is instructive to look at the interference effects between the 
various simultaneous and sequential contributions. To display these coher- 
ence effects, I choose to plot the scattering amplitude at zero degrees for the 
non-spin-fiip amplitude m p = m t = 1/2 (the only non-zero amplitude at 
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this angle). Fig. [3] plots all the simultaneous and sequential contributions 
from the different components listed in Table [TJ along with their coherent 
sums. We see that all the contributions to the simultaneous transfer are 
constructively coherent, as are all the contributions to the total sequential 
amplitude. This constructive coherence follows from the signs of the ampli- 
tudes in Table [TJ and reflects the significant pairing enhancement in 124 Sn. 
The total sequential and simultaneous amplitudes are not uniformly coher- 
ent with each other, however. This reflects the importance of the deuteron 
channel with its own specific optical potential. This is also shown by the 
rather different shapes of angular distributions in Fig. [TJ and the distinctive 
angle-dependent interference pattern. 
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